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Abstract 

In a class of the gauge-Higgs unification models the 4D neutral Higgs boson, 
which is a part of the extra-dimensional component of the gauge fields, becomes 
absolutely stable as a consequence of the gauge invariance and dynamically generated 
new parity, serving as a promising candidate for cold dark matter (CDM). We show 
that the observed relic abundance of cold dark matter is obtained in the 50(5) x U{1) 
model in the warped space with the Higgs mass around 70 GeV. The Higgs-nucleon 
scattering cross section is found to be close to the current CDMS H and XENONIO 
bounds in the direct detection of dark matter. 



What constitutes dark matter in the universe? How does the Higgs boson interact 

with other particles? [5] These are two of the most important issues in current physics. 
We would like to point out that these two issues are related to each other, and indeed that 
they reflect two sides of the fact that Higgs bosons can be stable and become dark matter. 

In the standard model (SM) of electroweak interactions, the electroweak (EW) symme- 
try is spontaneously broken by the Higgs field. However it is not clear at all if the Higgs 
boson appears as envisaged in the standard model. It is often argued that the naturalness 
and stability against radiative corrections to the mass of the Higgs field indicate the ex- 
istence of supersymmetry underlying the nature. However, there are other scenarios with 
the naturalness. The gauge-Higgs unification scenario is one of them, in which the 4D 
Higgs field, H{x), is a part of the extra-dimensional component of the gauge fields. It is 
identified as the fluctuation mode of the Aharonov-Bohm phase (Wilson line phase), 6h, 
in the extra-dimension. [6]- [8] 

It has been recently found in one of the gauge-Higgs unification models in the Randall- 
Sundrum warped space [9] that the Higgs mass ttih comes out around 50 GeV and the 
Higgs couplings with W and Z bosons, WWH and ZZH, and Yukawa couplings with 
fermions, Hipip, vanish at the value of 9h = i^^r where the effective potential for 9h is 
minimized. |10[ [TT] If this is the case, the Higgs boson would become stable. 

This prompts urgent questions. Does the Higgs couplings to the SM fermions and weak 
gauge bosons vanish to all order in perturbation theory? Is the value 6'// = ±|7r naturally 
achieved? Can the Higgs boson be absolutely stable? What physical consequences follow 
from stable Higgs bosons? 

We show in this paper that the WWH, ZZH, and Hifj if) couplings exactly vanish and 
the Higgs boson becomes stable, provided that certain conditions for the matter content 
are satisfied in the S0{5) x f/(l) gauge-Higgs unification in the Randall-Sundrum warped 
space. 

This leads to an important consequence. Higgs bosons become viable candidates for 
cold dark matter (CDM) in the universe. They are copiously produced in the very early 
universe. As the annihilation rate of Higgs bosons falls below the expansion rate of the 
universe, the annihilation processes get frozen and the remnant Higgs bosons become dark 
matter. 

The annihilation rates of Higgs bosons can be reliably estimated, once the masses of 
W, Z, and fermions are given. The precise value of the Higgs mass tuh, which arises as 
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quantum effects of Oh-, depends on the details of tlie model. In the model of ref. [10], mn 
is predicted to be ~ 50 GeV. It is shown below that the observed cold dark matter density 
is obtained with mu ~ 70 GeV in the 5*0(5) x f/(l) model. 

The distinctive feature in the gauge- Higgs unification is that the Higgs field H{x) corre- 
sponds to the fluctuation of the Aharonov-Bohm (AB) phase, Oh-, in the extra-dimensional 
space, which is given by the phase of the path-ordered Wilson line integral along a non- 
contractible loop C in the extra-dimension with a coordinate y, Pexp [iqa Jq dyAy[x, y)}. 
The effective Higgs interactions with the W , Z bosons and fermions in the standard model 
at low energies are summarized in the 50(5) x U{1) model as [11] 



a,b 
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Here k and L are two parameters specifying the Randall-Sundrum warped space. [9] The 
fundamental region in the fifth coordinate is < ?/ < L. The AdS curvature in the bulk 
region is given by A = — 6A;^. The warp factor = e^^ and the Kaluza-Klein (KK) 
mass scale ttikk = T^kzJ^^ are, typically, 10^^ and 1.5 TeV, respectively. The 5*0(5) gauge 
coupling qa is related to the 4D SU{2)l weak coupling g hy g = gA/ VL- 

The effective potential Vcs{Oh) is dynamically generated at the one-loop level, and 
is finite. [6] The mass functions 171^(9^), mziOn), and iti^ij^Oh) appear at the tree level, 
which include contributions from heavy KK excited states in intermediate states. [TT| [T2| [13] 
The value of the AB phase 6h is dynamically determined by the location of the global 
minimum of Vcs{6h)- There the fermion fields are diagonalized in the mass basis such that 

TT^ab = ^a^ab- 

In the 5*0(5) x U{\) model in the RS warped space [T^ [T2| [T5] the orbifold boundary 
conditions 





X, -y) = P, / (x, yj + y)p-' , (j =0,1), (2) 



are imposed for 5*0(5) gauge fields, where {yo,yi) = (0,/^). Here P/s are elements of 
50(5). With the parity matrices Pj = PJ'^'^ = diag (— 1, — 1, — 1, — 1, 1) in the vectorial 
representation, ^0(5) is broken to 50(4) ~ SU{2)l x SU{2)r. The 4D Higgs fields ^^(x) 
(a = 1, ■ ■ ■ ,4) appear as zero modes in the 5*0(5)/5*0(4) part of Ay{x,y). In components 
Ay{x,y) = (pa{x)ho{y) + ■■■ where /io(l/) is the zero mode wave function given below. 



An SO{A) vector (pa forms an SU{2)l doublet corresponding to the Higgs doublet in the 
standard model. 0i, 025 03 are absorbed by W and Z, whereas 04 corresponds to the 
physical neutral Higgs field. 

Let us denote the generators of S0{5)/ S0{4) by (a = 1, ■ ■ ■ ,4). In the vectorial 
representation (T4c)a6 = {i/ V^){Sa5Sb4 — Sa^Sb^), whereas in the spinorial representation 
T^^p = (l/2v^)/2 ® ri.[Il] The difference in normalization Tr T^T^^^ = 2TrT,^pT4 = 6ab 
becomes important in the subsequent discussions. 

The fifth dimension in the Randall- Sundrum warped space has topology of the orb- 
ifold S^/Z2. Its metric satisfies gMN{x,y + 2L) = gMN{x,y) = gMN{,x,—y). Conse- 
quently the AB phase Oh along the fifth dimension is given by exp{|6'//(2A/2T^)} = 
expjz^fA Jq dy{Ay)}. The relevant part of the gauge potential becomes 



Ay{x, y) = dH{x) ■ J ho{y) ■ + • • • , (3) 

where ho{y) = [2k/ {zl - l)Y/'^ e^^y {0 < y < L) and 9h{x) is defined in ([I]). ho{y) is 
normalized by dy e"'^'''^ ho{yy = 1, and is given, outside the fundamental region, by 
ho{-y) = ho{y) = ho{y + 2L). 

6h is a. phase variable. All the functions of 6h in the effective interaction in ([T]) are 
periodic with a period 27i. This periodicity follows from the large gauge invariance.|8], 
[T6l [T7] Given the boundary conditions ([2]), there remains the residual gauge invariance 
Am A\,j = QAm^~^ — {i/gA)^dM^~^ which preserves the boundary conditions. In 
general, new gauge potentials satisfy new boundary conditions 

X, y, -y)= P;A\,{x, y, + y)P;^' - — P; OmP;-' , (4) 

9a 

where = Q{x,yj — y)Pjfl{x,yj + y)~^. The residual gauge invariance is defined with 
fl{x,y) satisfying Pj = Pj. 

Consider a large gauge transformation Q}^'^^^{y) = exp |m dy \J Ak/[z\ — 1) /io(y) ■ 
T^} in the spinorial representation in which P^^ = /2 ® T3. It shifts 9h{x) to 0'u{x) = 
Oh{x) — ot- It is straightforward to see that the condition P'^^ = Pj^^ is satisfied if a = 27rn 
(n: an integer). In other words all physical quantities must be periodic in 6h with a period 
27r. 

Fermions in the vector representation of S0{5), for instance, obey \E'(x,yj — y) = 
Pj"^^l^^{^iyj + y)- Under a gauge transformation VL^^'^^°[y) with a = tt, one finds that 
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^p,vcc^ p,vcc^ = (Po™',Pi™') whereas (P^^^, Pf p) = (Pq'p, -Pi'p). Hence, if there are no 
fermions in the spinor representation of 5*0 (5), there appears enhanced symmetry. In this 
case all physical quantities become periodic in Oh with a reduced period vr. In the model 
of ref. [TD] bulk fermions appear only in the vector representation, thereby this condition 
being satisfied. Brane fermions located at y = are not affected by the transformation as 

There is mirror reflection symmetry in the extra dimension. The action in the 
RS warped spacetime is invariant under {x''^,y') = {x^,—y), AM{x,y) 

A'j^.j{x',y') = (Af,, -Ay){x,y), and ^{x,y) -> ^'{x',y') = ±Y^'^{x,y). The orbifold bound- 
ary conditions are preserved under this transformation. Since /io(— = ho{y), this implies 
that the theory at low energies is invariant under 9h{x) — > dni^') = —dnix)- 

In a class of the 5*0(5) x U (1) gauge-Higgs unification models in the warped space which 
contains fermions only in tensorial representations, but not in spinorial representations, of 
5*0(5), one can draw an important conclusion about the couplings of Higgs bosons. It 
follows from the enhanced gauge symmetry and mirror reflection symmetry that 

Kfr(^/f + tt) = Kff(^H) = Kfr(-^//) , 
^w,z{(^H + vr) = ml^zi^^u) = m'^w^zi-^n) , 

m^biOn + vr) = -m^M = m!^,{-9H) • (5) 

'^vk(O) = TTLziO) = as the EW symmetry is recovered at 6h = 0. The set of fermion 
masses {—m^f^} gives the same physics as the set {m^f^} does. The relative signs in the 
equalities for have been fixed by explicit evaluation at the tree level. [TT] 

It has been shown in ref. [10] that Kff(^H) is minimized precisely at 6h = i^vr due to 
the contribution from the top quark in the RS spacetime. The relations in ([5]), then, imply 
that all of the functions Vcs{Oh), t^w^z^^h), and rn^f^^On) satisfy a relation F{^'k + fj^^H) = 
P(|7r — f^^H). They are even functions of H when expanded around 9h = i^vr. 

It follows that all odd-power Higgs couplings H'^^+\ H^^+^WlW>', H^^+^Z^Z'', and 
H'^^^^ip a'ipb, vanish. In particular, the vanishing WWH and ZZH couplings signal signif- 
icant deviation from the standard model. Even if ttlh < m^y, the LEP2 bound for the 
Higgs mass is evaded. Further, Yukawa couplings for Hip aipb operators vanish identically, 
too. 

The derivative couplings such as 9^^h^7^7^^ and d^OnZ^^d'^Z'^ — d'^Z'^) are all for- 
bidden by the mirror reflection symmetry. We observe that the effective interactions at low 
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energies are invariant under H{x) —H{x) with all other fields kept intact at 9h = ±|7r. 
We call it the if-parity. Among low energy fields only the Higgs field is if-parity odd. The 
Higgs boson becomes stable, protected by the if-parity conservation. We stress that the 
if -parity has emerged dynamically, unlike in the models of refs. [T8| [T9] where an additional 
Higgs doublet with odd parity is introduced by hand. 

The mass functions are evaluated in the RS space. It is found in refs. |10 [ [TT | [T2] that, 
to a good approximation, 

mw0H) ~ cosdwrnziOn) ~ ^gfusmen , 

^a{(^H) ~ AaSin^H , (6) 

where 6*1^ is the weak mixing angle and the fermion mass matrix has been approximated 
by a diagonal one m^^^ = m^5ab- If a fermion belongs to spinor representation of 5*0(5), 
one would obtain mf ~ Aa sin \6h- As 6h = ^tt, one finds that my/ ~ \gfH and mf ~ Aq. 
The value of fn is given by fu ~ 246 GeV. We note that this differs from the vev of the 
Higgs field fuOH- 

Inserting ([6]) into ([I]), one finds the various Higgs couplings; 

~ -im^wK^'' + cos' ^ -Y, ai'a cos ^ . (7) 

The WWHH couphng is given by ^g^WjyV^H"^ , which is (—1) times the coupling in the 
standard modelll] This coupling includes contributions coming from tree diagrams contain- 
ing KK excited states Wn of W in the intermediate states with two vertices WWnH. [13] The 
ipipH'^ coupling is given by {ma/2f'fj)il) a'^^aH'^ ■ It is generated by two vertices ipipnH where 
ipn is the n-th KK excited state of ip. One comment is in order. The approximate formula 
for the fermion mass function in ([6]) may need corrections, depending on the details of the 
model. The symmetry property leads, in general, to m{6H) = X]^o^2n+i sin[(2n -|- 1)6'^^]. 
Accordingly the ipipH"^ coupling constant may be altered. 

Gauge-Higgs unification models under consideration are characterized with two param- 
eters fn and tuh at low energies. In a minimal model in the RS warped space, fn is fixed 
around 246 GeV by mw ~ \gfH- The value of m//, on the other hand, depends on the 
details of the matter content in the models. In the following numerical analysis, we fix 
fn = 246 GeV, whereas uih is treated as a free parameter. 
^We use diag.( — h ++) as 4D Minkowski metric. 
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With all the Higgs couplings at hand, one can estimate the annihilation rates of Higgs 
bosons in the early universe to determine its relic abundance as the cold dark matter. A 
rough estimate may be made with the following formula: VLfjh? ~ 3 x 10~^^ cm^s~^/(crf ), 
where VLh is the present mass density of the Higgs boson normalized by the critical density, 
h denotes the Hubble constant in units of 100 kms~^ Mpc~^, and {av) is the thermal average 
of the total annihilation cross section of the Higgs bosons multiplied by the relative velocity. 
The present mass density of cold dark matter is determined by WMAP collaboration as 
I^cdm/^^ = 0.1131 ± 0.0034.[2] 

Suppose that the Higgs mass is sufficiently smaller than mw- The dominant annihila- 
tion process is HH bb, and the abundance based on the above formula turns out much 
larger than the WMAP value. If the Higgs boson is heavier than W, HH W^W~ 
dominates, and the relic abundance turns out much smaller than the WMAP value. Thus, 
we expect that a Higgs boson slightly lighter than the W boson explains the cold dark 
matter abundance observed by WMAP. 

To determine a favored Higgs mass precisely, we have employed a more elaborated for- 
mula to evaluate the relic abundance [3],[20lll]. The annihilation rate per unit particle num- 
ber density av is expanded in a non-relativistic manner as av = a+bv'^+0{v'^), and the relic 
abundance is given by Qnh"^ — 2.82x 10® Yqo (mj^/GeV) where Y^^ = 0.264 gy^mpimH{a+ 
3{b — a/A)/xf}/xf. The freeze-out parameter Xf = niH/Tf, where Tf is the freeze-out tem- 
perature, is determined by Xf = ln[0.0382mpi(a + 66/x/')c(2-|-c)mi://((7*x/)"'^/^] with c ^ 0.5. 

We take the following annihilation modes into account: bb, V^*^V^*\ rf and cc. Here 
V denotes W or Z, and V* means a virtual gauge boson that eventually goes into a pair 
of a fermion and an anti-fermion. Since we are interested in the threshold region of the 
WW final state, the inclusion of 3-body and 4-body final states via virtual gauge boson(s) 
is mandatory. 

The cross section of the annihilation process into a final state X, a{HH X), is 
obtained from the decay rate into X, r{H X), in the standard model [3], by replacing the 
relevant vertices in the standard model by those in the present model, with an appropriate 
change in kinematical factors. For example, the cross section for HH — > // is given by 
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Higgs mass (GeV) 

Figure 1: Thermal relic density of Higgs boson DM with fn = 246 GeV. The solid (black) 
curve is obtained by the semi-analytic formulae, whereas the dashed (green) curve by the 
micrOMEGAs. The horizontal band is the WMAP data VlcDuh^ = 0.1131 ± 0.0034. 

where ^/s is the center-of-mass energy, /5j = a/1 — 4m^/s, and iV^ is the number of colors. 
The annihilation rate for HH y(*)v^(*) is given by 

(^VMVM = ^ dQiDviql) dql Dv{ql) , (9) 

Jo Jo 

where Dv{q^) = myTv/ {{q"^ —rriy)^ +myTy^ , and ctq is the annihilation cross section into a 
pair of intermediate virtual vector bosons, given by (Tq = {Gps/167rPi) 5yA^/^[A + 12g^g|/s^] 
with Sw(Z) = 2(1), and A(gf,g^; s) = s~^{{s - qf - g|)2 - Aqlq^}. 

Fig. [1] shows our numerical results for the relic abundance of the Higgs CDM. The 
solid (black) curve is obtained by the semi-analytic formulae with all the aforementioned 
channels included in the total annihilation rate. The horizontal (blue) lines near 0.1 indicate 
the allowed range by WMAP. The favored Higgs mass is found around 70 GeV in a rather 
narrow range. The freeze-out parameter is x/ ~ 21 in the favored region, which corresponds 
to T/ ~ 3 GeV. 

At niH = 70 GeV, the annhilation rates of the relevant modes in the non-relativistic 
limit are av = (7.3, 11, 1.5) x lO'^^cm^/s for bb, W^*^W^*\ and Z^*^Z^*^ respectively. The 
mode is larger than the bb mode even below the WW threshold, which confirms 
the importance of the iy(*)l4^(*) mode. All the remaining modes including the gg mode 
have smaller rates. 

Importance of the 3-body and 4-body final states through HH below the 

threshold of VV pair can be seen by comparing the solid (black) curve with the dashed 
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(green) curve, which was obtained using the micrOMEGAs 2.2 |21j with two-body final 
states only. The WW channel opens even below the threshold due to the thermal energy. 
The micrOMEGAs including only 2-body final states gives mn — 75 GeV as a favored Higgs 
mass. The relative contributions of the bb and WW modes are 34% and 61%, respectively. 
The result obtained from the semi-analytic formulae agrees with the result from the mi- 
crOMEGAs with only 2-body final states well below the WW threshold. However, near the 
WW threshold, the relic density from micrOMEGAs with 2-body final states in the HH 
annihilation yields substantially larger Qnh"^ than the semi-analytic treatment, although 
it includes thermally allowed 2-body final state HH WW. Near the WW threshold, it 
is important to include the virtual W effect correctly in order to get an accurate behavior 
of the relic density across the threshold. 

If Higgs bosons constitute the cold dark matter of the universe, they can be detected 
by observing Higgs-nucleon elastic scattering process, HN — > HN. The relevant part of 
the effective interaction ([7]) is CcH = {H"^ I'^fh) Ylf'^ff f- 

To evaluate the direct detection rate one needs to incorporate QCD corrections [22] . 
After integrating out the heavy quarks, we obtain the effective Lagrangian at a hadronic 
scale: 

7-7-0 

Air' 



r ^ 



(10) 



_q=u,d,s 

where c, b and t quarks are integrated out, and denotes the gluon field strength. 
This effective Lagrangian leads to the following effective Higgs-nucleon coupling: 

C„. ^ H H-NN , (11) 

where h = E,.„,<i,. o-'^i Wm,qqm = m^/,". The relation {JV|(cv./8ir)GG|]V> = 
— (mAr/9)(l — /at) has been used. With this coupling, the spin-independent (SI) Higgs- 
nucleon scattering cross section is evaluated to be 

1 ^ 2 + 7/^ V m% 



in the non-relativistic limit. There is a considerable uncertainty in the value of Jn stemming 
from that in the ct^at term. The value /at quoted in ref.'s and [23] are /at ~ (0.31 ~ 0.41) 
and /at ~ (0.20 ~ 0.45), respectively. On the other hand, the recent lattice calculation 
gives a smaller value for and /at ~ 0.07 [24j. We choose /at = 0.3 and 0.1 for the 
purpose of illustration in the following. 
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Figure 2: Spin-independent HN scattering cross sections as functions of Higgs mass along 
with the bounds of CDMS II (dashed, green) and XENONIO (dasd-dotted, blue). The 
solid (black) and the dotted (red) curves are for two different values of /at = 0.1 and 0.3, 
respectively. 



In fig. [21 we show the Higgs-nucleon scattering cross sections as functions of Higgs mass 
for two different values of /at = 0.3 and 0.1 with fn = 246 GeV. The present experimental 
upper bounds for the spin-independent WIMP-nucleon cross sections from CDMS II [25] 
and XENONIO [26] are depicted as well. For mu = 70 GeV, the experimental bound is 
c^si ^ 5 X lO'^^'cm^ at 90 % CL, whereas our prediction is asi — (1.2 — 2.7) x lO^^'^cm^ 
for = (0.1-0.3). 

This does not necessarily mean that the present model is excluded. The direct de- 
tection signals are proportional to the local density of cold dark matter, po? which has 
not been measured by experiments. In most cases including CDMS II and XENONIO, 
the experimental bounds are derived under assuming po = 0.3 GeV/cm^. For a spherical 
and smooth halo, po = 0.2 — 0.6 GeV/cm^ seems a reasonable range [3l H]. Taking the 
lower value relaxes the constraint. Further it has been argued that po can be as small as 
0.04 GeV/cm^ for non-smooth distribution of dark matter in the Galactic halo [2^, which 
makes mn = 70 GeV consistent with the data. On the theoretical side, as mentioned below 
eq. ([7]), the HHff coupling constants may be reduced if the mass function mf{6H) has 
more general ^^-dependence. In this case the Higgs-nucleon cross section is decreased, 
while the Higgs relic abundance, which depends on both the HHhh and HHWW cou- 
plings, is kept unchanged by appropriately increasing mu- To summarize, it is premature 
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to exclude the Higgs dark matter scenario based on the current bounds from CDMS II and 
XENONIO. 

Signals from pair annihilation of Higgs bosons in the Galactic halo into two 7's or a 
7 and a Z boson may be seen, too. Expected signals are two (nearly) monochromatic 
gamma lines of = mni— 70GeV) and -E-y = mn — m'^/{4'mH){— 40GeV). Their rates 
are estimated to be cr^^w|„^o — 4.3 x 10~^^cm^/s, and cryzv\v~,o — 5.4 x 10~^^cm^/s 
for rriH = 70 GeV and fn = 246 GeV. Comparing these rates with those in the inert 
doublet model [19], we expect that such monochromatic gamma rays could be observed by 
FERMI Gamma-ray Space Telescope (formerly GLAST). However the signals will be less 
pronounced in our model compared with the inert doublet model, since the annihilation 
cross section is dominated by HH bb,W^*^W^*^ in our case, which yield secondary 
photons with a continuum spectrum. 

In this letter we have shown that the Higgs boson becomes absolutely stable in a class 
of the gauge- Higgs unification models. The stability of the Higgs boson is protected by a 
new dynamically emerging parity, if-parity. Higgs bosons become the cold dark matter in 
the universe. The observed dark matter density is obtained with mn ~ 70 GeV. Although 
the direct detection rate for the Higgs-nucleon elastic scattering is found slightly above 
the current upper bounds, the prediction in the gauge-Higgs unification model may be 
consistent in view of many uncertainties involved. If the Higgs boson is absolutely stable, 
the way of finding Higgs bosons in collider experiments must be scrutinized. Higgs bosons 
appear as missing energies and momenta. We shall come back to these points in future. 
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